An Alternative Parametrization of the Pion Form Factor 
and the Mass and Width of p{770)* 



Peter Lichard^'^ and Martin Vojik^ 

'institute of Physics, Silesian University in Opava, Bezrucovo nam. 13, 74-6 01 Opava, Czech Republic 
Institute of Experimental and Applied Physics, Czech Technical University, Horskd 3/a, 120 00 Prague, Czech Republic 

(Dated: February 2, 2008) 

In order to reveal possible mass shifts of the vector mesons in a dense strongly interacting medium 
presumably created in high energy heavy ion collisions it is necessary to know their free masses 
reliably. The p(770) mass quoted in the last two editions of the Review of Particle Physics is 
significantly larger than the values quoted in previous editions. The new value is mostly influenced 
by the results of recent experiments CMD-2 and SND at the VEPP-2M e'''e~collider at Novosibirsk. 
We show that the values of the mass and width of the p(770) meson measured in the e+e~ tt+tt" 
annihilation depend crucially on the parametrization used for the pion form factor. We propose 
a parametrization of the p(770) contribution to the pion form factor based on the running mass 
calculated from a single-subtracted dispersion relation and compare it with the parametrization 
based on the formula of Gounaris and Sakurai used recently by the CMD-2 collaboration. We show 
that our parametrization gives equally good or better fits when applied to the data of the CMD- 
2, SND, and KLOE collaborations, but yields much smaller values of the p(770) mass, consistent 
with the photoproduction and hadronic reactions results. Our fit to the KLOE data becomes 
exceptionally excellent (confidence level of 99.88%) if an energy shift of about 2 MeV in the p-ui 
region is allowed. 



I. INTRODUCTION 

The vector meson p(770) plays an important role 
in many phenomena of particle and nuclear physics. 
Possible modification of its mass and width in dense 
strongly interacting systems presumably created in high 
energy heavy ion collisions is now widely discussed [1- 
11]. To study this effect quantitatively, reliable informa- 
tion about the "vacuum" values of these parameters is 
desirable. 

Inspecting recent editions of the Review of Particle 
Physics [12-19] we find: (i) whereas the values of /o(770) 
mass quoted until 2002 are compatible within experimen- 



TABLE I: Mass and width of p(770) according to PDG 



year 


mp 


Pp 




(MeV) 


(MeV) 


1992 [12] 


768.1 ±0.5 


151.5 ± 1.2 


1994 [13] 


769.9 ±0.8 


151.2 ± 1.2 


1996 [14] 


768.5 ±0.6 


150.7 ± 1.2 


1998 [15] 


770.0 ±0.8 


150.7 ± 1.1 


2000 [16] 


769.3 ±0.8 


150.2 ±0.8 


2002 [17] 


771.1 ±0.9 


149.2 ±0.7 


2004 [18] 


775.8 ±0.5 


150.3 ±1.6 


2006 [19] 


775.5 ±0.4 


149.4 ± 1.0 



*This paper is dedicated to the late Julia Thompson, who drew the 
attention of one of us (P. L.) to the experimental program of the 
Budker Institute of Nuclear Physics at Novosibirsk. 



TABLE II: Mass of p measured in different processes, selected 
and averaged by PDG in 2006 [19] 







mp± 




(MeV) 


(MeV) 




775.5 ± 0.4 


775.4 ± 0.4" 


Photoproduction 


768.5 ± 1.1 


/ 


Hadronic reactions 


769.0 ± 0.9 


766.5 ± 1.1 



'Includes also r decay results 



tal errors, an abrupt rise appears in 2004 (see Table I); 
(ii) the p° mass average calculated from the e+e~ an- 
nihilation experiments is much higher than the averages 
calculated from hadronic reactions and photoproduction. 
The figures taken from the 2006 Particle Date Group 
(PDG) tables are shown in Table II. 

The discrepancy between the e^e" annihilation mea- 
surements of the p° mass and the results of photoproduc- 
tion and hadronic experiments is very troublesome. We 
know of no single theoretical prediction or explanation of 
such a phenomenon. 

In the literature one can find many articles that dealt 
with the determination of the p(770) parameters from 
the experimental data on the pion form factor [20-30]. 
Their results on the p'^ mass are smaller, in some cases 
substantially, than the value preferred now by the PDG. 
It is not clear what has led the PDG to not using them 
"for averages, fits, limits, etc." [19]. The PDG aver- 
age builds on five experiments, three of which explored 
the e'^e~ annihilation into a pair of charged pions. Two 
of those three experiments utilized the Gounaris-Sakurai 
(GS) parametrization [31] when fitting the data. 

As we discuss in detail in the next section, the valid- 
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ity of the GS formula is limited because it ignores two 
important effects. Namely, the opening of additional de- 
cay channels at higher energies and the structure effects 
in the hadronic vertices manifesting themselves as the 
momentum dependent strong form factors. These effects 
directly influence the imaginary part in the denominator 
of the GS formula. This itself would not have serious im- 
plications for the analysis of the form factor data, which 
are usually restricted to an narrow interval around the 
nominal p mass. But because of the close connection be- 
tween the denominator's real and imaginary part implied 
by the analyticity, the real part is also influenced, even 
at the energies around the nominal p mass, which are 
important for data fitting. 

Already in the pioneer era of the p-meson physics, 
when the p mass and width were determined mainly 
by analyzing the dipion mass distribution in the re- 
action Tip ttttN, Pislit and Roos [32] showed that 
the discrepancies among the results of different experi- 
mental groups are caused by different parametrizations 
these groups used. They also discussed the discrep- 
ancy between the hadronic measurements and the re- 
sults of the first e+e~ annihilation into two pions ex- 
periment [33]. From that time the parametrization de- 
pendence, mainly in connection with the pion form factor 
parametrization, has been discussed in many papers, see, 
e.g., Refs. [24, 26, 28, 30, 34, 35]. 

The aim of this paper is to show that the uncertainties 
connected with the choice of the form factor parametriza- 
tion induce much larger errors of the p resonance param- 
eters than the statistical and systematic errors of the 
data. For that purpose we fit three different data sets 
with two different parametrizations and compare the re- 
sulting masses and width of the meson. One of the 
parametrizations is that used recently by the GMD-2 
collaboration [36, 37]. Its principal element is the GS 
formula. The other parametrization [38] is based on a 
single-subtraction dispersion relation fed by the imag- 
inary part taking into account several decay channels 
which open when the p-resonance goes "off-mass-shell" 
towards higher masses. All decay widths are calculated 
taking into account the strong form factors. 

The paper is organized as follows. In Sec. II we de- 
scribe the two versions of the p'^ resonance term in the 
pion form factor. The parametrization of the sublead- 
ing contributions coming from the p(1450), aj(782), and 
(/ii(1020) is presented in Sec. III. The results of fitting the 
CMD-2 [37], SND [39], and KLOE [40] data sets by the 
two parametrizations are shown and discussed in Sec. IV. 
Section V contains summary and concluding remarks. 



in the following form 



M2(0) 



M^{s)-s-impTp{s)' 



(1) 



where the running mass squared M'^[s) is a real function 
satisfying the conditions 



and 



ds 



= 0, 



(2) 



(3) 



rUp is the nominal mass of the p meson, and the real func- 
tion Tp{s) is the total decay width of a p° meson with 
mass ^/s. It follows from microcausality that the pion 
form factor F{s) is an analytic function in the whole com- 
plex s-plane with a cut running from s = Arri^ to infinity, 
satisfying the condition F{s*) = F*{s), and growing no 
faster than a finite power of s for |s| oo (in QCD, |-F(s)| 
even decreases [41]). These properties of F{s) imply the 
analytic properties of the denominator in Eq. (1). 



A. Gounaris-Sakurai formula 

Gounaris and Sakurai [31] assumed a low-energy 
parametrization of the phase shift in the P\\ wave of 
the TTTT elastic scattering inspired by work of Chew and 
Mandelstam [42] and used the relation between the pion 
form factor and the tttt elastic scattering amplitude in 
the Pii channel [43]. After defining the resonance pa- 
rameters in terms of the Pn phase shift they obtained a 
parametrization of the pion form factor. Below, we write 
it in a little different form, adopting partly the notation 
of Ref. [36] and referring to our Eq. (1). Firstly, on the 
upper branch of the cut we have 



(4) 



where Fp — V p{m?p) is the nominal width of the p(770) 
meson, whereas Fp(s) = for s < 4m^. Secondly, 



Ml{8)=ml + f{s), 



(5) 



II. TWO PARAMETRIZATIONS OF THE p(770) 
CONTRIBUTION TO THE PION FORM FACTOR 

We write the contribution of the intermediate state 
with the p^ meson to the pion electromagnetic form factor 



where 

2r m? f 

+ {ml-Aml)h'{m'^{ml-s)], (6) 
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TT V « 2to„ ' — 



/i(s) = < 



< s < 4ml, 



1 — 4 arctan 



/ 4m J- 



s = 0, 



s-4mj i„ ^/^+■^4m^-s 



To get more insight into the physics behind the GS 
parametrization let us utilize that M^(s) and Tp{s), 
which arc related to the p-mcson contribution to the pion 
form factor by Erj. (1), should be the bomidary values of 
the real and imaginary part, respectively, of a function 
analytic in the cut s-plane. It is easy to check that the 
quantities (4) and (5) are related through the double- 
subtracted dispersion relation 



M'^{s) = Ao + Ais ■ 



OO 

TT J S'^{s' 



s) 



(7) 



The appearance of the two subtraction constants Aq and 
Ai guarantees that the conditions (2) and (3) can be 
satisfied. Of course, the relation (7) is not specific for the 
GS parametrization, but is more general. It enables us 
to find the running mass squared M^(s) corresponding 
to any reasonable (i.e., not growing too fast) choice of 
the total width rp(,s). It is important to realize that 
the values of Mp{s) at any s, in particular those for low 
s < 1 GcV^ , which are used in the form factor fitting, are 
given by the values of rp(s) at all s. To get reliable Mp{s) 
it is therefore essential to have rp(s) that is physically 
sound not only in the s-region vital for the fitting, but 
also beyond it. 

The GS parametrization of the p{770) contribution 
to the pion electromagnetic form factor ignores possi- 
ble structure effects in the p^tt+tt" vertex (the strong 
form factor). They were taken into account by Vaughn 
and Wali (VW) [44] in the form of a phenomenological 
cutoff function. When that function is set to unity, the 
VW parametrization reduces to the GS one. The VW 
approach has not become as popular as that of GS, al- 
though it is closer to reality, because the outcome of dis- 
persion relation (7) cannot be expressed in an analytic 
form, similar to (6). 

A different attempt to bring the GS-type formulas 
closer to reality was done by the authors of Rcf. [30]. 
They included also the contributions from the K~^K~ 
and K^K'^ states to the imaginary part of the denom- 
inator of the pion form factor (1) and obtained its real 
part in an analytic form using the dispersion relation. No 
strong form factor was included. Also this improvement 
of the GS formula passed unnoticed by the experimen- 
talists fitting their data. 



We should also note that the original derivation of 
the GS formula brings a one-to-one correspondence be- 
tween the Pii scattering length and effective radius on 
one side, and the p(770) resonance position and width 
on the other. Using the GS formula for the contribution 
of higher resonances [p(1450), p(1700)] to the pion form 
factor, as in [36, 37, 45, 46], thus contradicts the spirit of 
the GS derivation as it implies incorrect values of the Pn 
scattering length and effective radius. In the language of 
dispersion relation (7) it corresponds to the assumption 
that the two-pion mode is the only decay mode at any 
y/s. This is imrcalistic because already for the p(1450) 
with nominal mass there are several other modes allowed 
by the energy-momentum, angular momentum, and in- 
ternal quantum numbers conservation. 



B. Our running mass approach 

We build our alternative parametrization on the recent 
work by one of us [38], to which the reader is referred for 
technical details. The numerical values of the constants 
used in [38] arc updated using the 2002 edition of the 
Review of Particle Physics [17]. 

There are four items in which our parametrization dif- 
fers from the GS one: 

(i) We consider several contributions to the total decay 
width of the p meson, which become operational when its 
mass ^/s exceeds respective thresholds 



rp(s) 



- /9— »7r7r 
■ p^KK 



/9— s-r/TTTT 



is). 



(8) 



(ii) The partial decay width of the two-pion decay, 
which provides the essential contribution to (8) is cal- 
culated from the standard interaction Lagrangian con- 
taining the first derivative of the pion field. The result 



IS 



r(s) = r, 



m„ 



Am. 



3/2 



mt 



Ami J 



(9) 



(iii) The strong form factors are introduced into all 
purely hadronic vertices. They were taken from the fiux- 
tubc-brcaking model of Kokoski and Isgur [47]. This 
leads to the modification of the partial decay widths com- 
posing (8) and calculated originally as if the participating 
mesons were elementary quanta. For example, the two- 
pion width (9) is modified by the following multiplication 
factor 



exp • 



24/32 



where the value /? = 0.4 was taken from the Kokoski and 
Isgur [47] analysis. 

(iv) As now the the total width (8) docs not rise so fast 
as in the GS approach, we may write a once-subtracted 
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Gounaris-Sakurai 



■Gounaris-Sakurai 



- Running mass [20] 



-Running mass [30] 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

s (GeV2) 




FIG. 1: Comparison of tlic running mass squared of p- meson 
in the Gounaris-Sakurai and our approach 



FIG. 2: Comparison of the energy dependent width of p- 
meson in the Gounaris-Sakurai and our approach 



dispersion relation 

oo 

^ TT J S'{s' — S) 



(10) 



4mJ 



The subtraction constant Bq is chosen to ensure that 
the condition Mp(rn^p) = nip (2) is satisfied. Now it 
is not possible to satisfy the condition on derivative (3) 
automatically. We consider this an advantage, because 
Eq. (3) can thus serve as a test that the energy dependent 
total decay width (8) was chosen correctly. For details, 
see Ref. [38]. Let us note that our choice of the sub- 
traction point s = differs from s = m^, used in Refs. 
[38, 48] . The small disadvantage lying in the necessity to 
evaluate the subtraction constant Bq is generously com- 
pensated by much faster numerical evaluation of the dis- 
persion integral. In the following, we are going to vary 
nip and Tp in a search for the best fit to the experimental 
data. It means that we will have to calculate the running 
mass (10) in the experimental points again and again. 
The computation speed is therefore a decisive factor. 

The comparison of the basic quantities of the GS and 
our approach is shown in Fig. 1 (running mass) and Fig. 2 
(total width). The values nip = 771.1 MeV/c^ and Tp = 
149.2 MeV [17] were chosen for this purpose. 



III. OTHER CONTRIBUTIONS TO THE PION 
FORM FACTOR 

The parametrization of the pion electromagnetic form 
factor should contain, in addition to the dominant contri- 



bution from the /9(770)*^, also the contribution from other 
resonances. In the ciiiergy range of the e+e^ tt+tt^ ex- 
periments we are interested in, they include p' = p(1450), 
w(782), and perhaps also ^(1020). We write the pion 
form factor in the form 

F{s) = \Fp{s)[l + 5A^F^{s)\+l3Fp.{s)+^F4,{s)] 



X (1 + /3 + 7)- 



(11) 



where not only F{s) but also the individual contributions 
to it are normalized to unity at s = 0. Parameters 0, 7, 
and S are complex and should be determined, together 
with nip and Tp by minimizing constructed from the 
data and parametrizations in the standard way. Equa- 
tion (11) is inspired by Eq. (8) from [36] and reduces to 
it for 7 = 0. We will use it in conjunction with both 
GS and our running mass parametrizations. In one case 
(CMD-2 data) allowing non-vanishing 7 did not bring a 
considerable drop in x^- But for the other two data sets 
(SND, KLOE) the cj) contribution improved the fits. 

The CMD-2 collaboration [36, 37] used the GS 
parametrization also for the p(1450) contribution and we 
will do the same when applying their formalism. But 
when using our parametrization, we utilize the fixed-mass 
formula 



m 



p' 



m^p, - s - imptTpt{s) 



(12) 



with several options for the total width of p': 

1. Tp:{s) = rp(i45o); 
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2. Eq. (4) adapted for p(1450); 

3. r^'(s) = 0; 

4. Eq. (9) adapted for p(1450). 

The options 2 and 4 are equivalent to the assumption that 
the only decay mode of p(1450) is the dipion one. The 
option 3 mimics the situation when p(1450) disintegrates 
only through the modes with higher thresholds, the de- 
cay widths of which become negligible below s = 1 GeV^ 
mainly due to the shrinking phase-space. This is in con- 
formity with what was known already a long time ago, 
namely that the p' resonance is clearly visible in the e"*"e~ 
annihilation into four pions [49] whereas its appearance 
in the two-pion channel is less pronounced, see, e.g., Fig. 
2 in [46]. 

The results of the fits (x^, TOp(77o), ^ pijiQ)) showed 
only marginal differences among various options, with 
option 3 giving slightly better than the other options. 
The results of our approach presented below were ob- 
tained using just that option. 

As to the contribution of the a;(782) to the pion elec- 
tromagnetic form factor, we use 



mr. 



mf. 



im. 



(13) 



The following options are considered: 

1. Tt^{s) = r(^(782) above the s = 9m^ threshold; 

2. r^(s) = r„^3.(s); 

3. same as 2, but including strong form factors; 

4. r^{s) = Tu^^snis) + Ta,^^j{s) + r„^^+„- (s); 

5. same as 4, but including strong form factors. 

Following [36], the option 1 was used in conjunction with 
the GS parametrization of the p(770) and p(1450) con- 
tributions. The other options were explored as a part of 
our approach, but did not exhibit notable differences in 
the quality of the fits. 

In the case of the SND data, the contribution of the p-(f) 
interference is present in a statistically significant way, as 
judged from the decrease, with both parametrizations 
(GS, ours) which we explore. For data of the KLOE 
collaboration, the situation is more complicated. The 
details are reported in Sec. IV C. Our parametrization of 
the ^(1020) contribution to the pion form factor is 



F^is) 



ml 



ml - s- im^T^{s)' 



(14) 



We consider the following options: 

1. r0(s) = r0(io2o) above the s = threshold; 

2. T^{s)=T^^^-j^{sy, 

3. same as 2, but including strong form factor; 



4. r^(s) = T^^j^j^{s) + r^^3^(s) + r^^^^(s); 

5. same as 4, but including strong form factors. 

Given the narrow width of the (j) resonance, the high 

thresholds of its important decay modes, and the dis- 
tance of the experimental energies from the (p mass, the 
refinements 2 5 did not c;hange the results noticeably in 
comparison with the simplest option 1. 

IV. RESULTS 
A. Using the CMD-2 experiment data 

In 2002, the CMD-2 Collaboration published [36] their 
results of the precise measurement of the e+e^ — > tt+tt" 
cross section for 43 cm. energies and gave also the values 
of the pion electromagnetic form factor determined from 
them. After an error in the computer code was found, 
they published the updated values of the pion form factor 
[37], which we are going to explore here. The original 
paper [36] remains important as it contains details about 
the GS based parametrization, which the collaboration 
used in fitting the form factor data. 

We first tested our computer code using the uncor- 
rected CMD-2 data because in Ref. [36] a more complete 
information is provided than in the later paper [37]. Our 
results agreed with the values in Table 3 [36]. Then we 
used the corrected pion form factor data [37] and got the 
results shown in Tab. III. The mass and width of the 
p meson we obtained differ only very slightly from the 
CMD-2 values, differences are much smaller than the er- 
rors quoted in [37]. Our errors shown in Tab. Ill and 
elsewhere are only the statistical ones in the MINUIT 
[50] sense. 

Our running mass parametrization gives a little better 
fit to the CMD-2 data than the GS one. The results are 
shown in the last column of Tab. III. The p(770) mass 
and width we got are much lower that the corresponding 
quantities obtained by the GS parametrization. 

It is also interesting to compare the mean square radius 
of the charged pion, which is related to the derivative of 
the form factor at s = [51], with the experimental value 
< rl >=(0.439±0.008) fm^ [52]. 

The graphical comparison of our parametrization with 
the CMD-2 data is shown in Fig. 3. 



B. Using the SND experiment data 



In 2005, the results were reported [53] of the e+e~ — > 
TT+TT" cross section measurement with the spherical neu- 
tral detector (SND) at the VEPP-2M collider in Novosi- 
birsk. Later on, an error in software was found and the 
data corrected [39]. We use the data on pion form factor 
from their Table I. The experimental value of the charged 
pion mean square radius [52] is added to the fitted data 



6 



TABLE III: Mass and width of p(770) and the mean-square 
radius (msr) of the pion from fitting the CMD-2 data. The 
experimental value of the msr is 0.439±0.008 fm^ [52]. 





GS 


ours 


X^/NDF 


35.16/37 


32.18/37 


rup (MeV/c^) 


775.3 ±1.1 


767.08 ± 0.83 


Tp (McV) 


143.2 ± 2.5 


136.1 ± 2.2 


< ri > (fm^) 


0.416 ± 0.013 


0.434 ± 0.012 



50 




Vs (MeV) 



FIG. 3: Comparison of our parametrization of the pion form 
factor with the CMD-2 data [37] 

set. The comparison of the GS and our parametrization 
(both include also the contribution from the (f) meson) is 
reported on in Tabic IV. Wc also reproduce there the 
results of the updated [39] fit by the SND collaboration. 
The formulas used in their fitting are described in [53] 
together with three choices of the p' and p" parameters. 
Unfortunately, in the updated paper the provided infor- 
mation is scarce. It is not said which choice was used 
and what was the of the fit. Nevertheless, there is a 
good compatibility in both rUp and Tp between the GS 
parametrization and that used by the SND collaboration. 

The GS and our parametrization provide equally good 
fits, but the mass and width resulting from the latter are 
lower than and incompatible with those from the former. 

C. Using the KLOE experiment data 

The data on the cross section a{e~^e~ — > 7r+7r~7) at 
the collision energy W =1.02 GeV measured with the 
KLOE detector at the electron-positron collider DA$NE 



TABLE IV: Mass and width of p(770) from fitting the SND 
data. 

GS ours SND fit [39] 

X^/NDF 39.94/38 39.80/38 not given 

rUp (MeV/c^) 774.60 ± 0.52 764.23 ±0.43 774.6 ± 0.4 ± 0.5 
(MeV) 147.3 ± 1.0 140.29 ± 0.85 146.1 ± 0.8 ± 1.5 



50 




Vs (MeV) 



FIG. 4: Comparison of our parametrization of the pion form 
factor with the SND data [39] 



in Frascati were used in Ref. [40] to extract the square 
of the charged pion form factor |^V(s)|^ by the radiative 
return method [54]. Here, s < W"^ denotes the square 
of the invariant mass of the tt+tt" system. No attempt 
was made by the authors of [40] to fit their form factor 
in order to get the p meson parameters. We have per- 
formed the fit using both GS and our parametrizations. 
The results arc summarized in Table V. The of both 
fits is poor. Nevertheless, the of our parametriza- 
tion is about half that of the GS parametrization. An 
interesting observation is that not only OTir parametriza- 
tion, but also the GS one gives a p(770) mass substan- 
tially smaller than both Novosibirsk experiments. As be- 
fore, our parametrization yields smaller values of both 
the mass and width than GS. An attempt to improve 
the quality of the fits by including the p-(p interference 
somewhat lowered the -y^ for both parametrizations, but 
not enough to get the confidence level into an acceptable 
range. 

Figure 5 shows that the worst discrepancy of our 
parametrization with the data is located in the p-uj re- 
gion. The abrupt fall, which is present in Figs. 3 and 4, 
is missing. A more detailed analysis of the latter figures 
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TABLE V: Mass and width of p(770) from fitting the KLOE 
data. 





GS 


ours 


xVndf 


367.8/54 


192.6/54 


rup (MeV/c^) 


769.22 ± 0.32 


761.39 ± 0.25 


(MeV) 


144.71 ± 0.56 


140.07 ± 0.66 




500 600 700 800 900 1000 

Vs (MeV) 

FIG. 5: Comparison of our parametrization of the pion form 
factor with the KLOE data [40]. 



showed that the steep descent is caused by the p-u inter- 
ference term, which is very sensitive to the exact position 
of the (very narrow) cu resonance. The absence of the in- 
terference term in Fig. 5 suggests that the location of the 
to resonance in the KLOE data is different from the nom- 
inal one, which we assumed when doing the fit with both 
GS and our parametrization. It would mean that the en- 
ergy scale in the KLOE data contains a systematic shift, 
at least in the p-ui region. The KLOE collaboration uti- 
lized the radiative return method [54]. This means that 
the invariant mass of the 7r^7r~7 system (denoted as W 
in [40]) is fixed by the energy of the colliding e+ and 
e~ . But the invariant mass ^/s of the 7r"'"7r~ subsystem 
has to be calculated from the measured momenta of the 
outgoing particles and is therefore known less precisely. 

To scrutinize the possibility of a shift in experimen- 
tal values of y/s, we repeated the fits allowing the co 
mass also to vary. Both the GS and our parametriza- 
tion indicate that the best fit is obtained if the lu mass is 
smaller than the nominal one by about 2 MeV, see Ta- 
ble VI. The drop of in both parametrizations justifies 
including the lu mass among free parameters, compare 
Table VI with Table V. But the quality of the fit is very 




\/s (MeV) 

FIG. 6: Comparison of our parametrization of the pion form 
factor with the KLOE data [40] after decreasing the a;(782) 
mass parameter by 1.93 MeV. 



TABLE VI: Fitting the KLOE data, the a; (782) mass allowed 
to vary. 

GS ours 



X^/NDF 

(MeV/c^) 

nip (MeV/c^) 
rUp (MeV/c^)' 
Tp (MeV) 



202.6/53 
780.52 ± 0.33 

2.13 ±0.35 
769.45 ± 0.25 
771.58 ± 0.60 
146.39±0.51 



29.58/53 
780.72 ± 0.15 

1.93 ±0.19 
761.51 ±0.14 
763.44 ± 0.33 
142.56 ± 0.41 



"Difference from m„ = (782.65 ± 0.12) MeV [19] 

'The value resulting from the fit was increased by Am^ 



different for the two parametrizations. While the value of 
= 202.6 for the GS parametrization is still very high 
(almost four times the number of degrees of freedom) , our 
parametrization yields = 29.6, which means, consider- 
ing NDF = 53; the confidence level of 0.9962. Such per- 
fect fits are rarely encountered in the high energy physics. 
The improvement of the fit is clearly seen by naked eye if 
one compares the new parametrization (Fig. 6) with the 
old one (Fig. 5). 

Assuming that our analysis really signifies that a sys- 
tematic shift of is present in the data, we used the 
difference Am^, between the nominal mass of the a;(782) 
and that resulting from our fit to get the "corrected" 
p-meson mass, see Table VI. 

Unfortunately, our interpretation of the mass shift in 
the p-Lu mass region as an experimental effect is difficult 
to be explained in terms of a genuine miscalibration of 
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TABLE VII: As Tab. VI, but including the p-0 interference. 

GS ours 



X^/NDF 

(MeV/c^) 
Amu. (MeV/c^)" 
rup (MeV/c^) 
mp {MeV/d'f 
Tp (MeV) 



58.35/51 
780.15 ±0.17 

2.50 ±0.21 
770.29 ±0.16 

772.79 ± 0.35 
151.53±0.55 



25.67/51 
780.71 ±0.15 
1.94 ±0.19 
761.54 ±0.16 
763.48 ± 0.33 
142.78 ± 0.46 



"Difference from = (782.65 ± 0.12) MeV [19] 

''The value resulting from the fit was increased by Am^ 



TABLE VIII: Confidence level of the GS and our fits to vari- 
ous data sets 



GS 



ours 



CMD-2 

SND 

KLOE 

KLOE" 

KLOE* 



0.5555 
0.3840 
< 10-30 
2x10"^^ 
0.2234 



0.6943 
0.3899 
2x10"" 
0.9962 
0.9988 



"With the to mass shifted 
''With the uj mass shifted and • 



contribution included 



the KLOE drift chamber [55]. 

In an effort to reveal the origin of the mysterious mass 
shift wc considered the possibility that it was an ar- 
tifact of a wrong parametrization. Wo combined our 
parametrization of the contribution (1) to the pion 
form factor, which is based on Eq. (10), with various 
options for other contributions, see Sec. III. The ef- 
fect did not disappear and the value of the uj mass shift 
was very stable against all changes. Moreover, we have 
found that the inclusion of the (?!)(1020) contribution fur- 
ther improves the quality of the fit not only for our 
parametrization (where the quality was already excel- 
lent), but also very significantly for the GS parametriza- 
tion, see Tab. VII. 



V. SUMMARY AND COMMENTS 

The main results of this work are summarized in Tables 
VIII, IX, and X. 

Table VIII shows how well are the particular data 
sets (CMD-2, SND, KLOE) fitted by two parametriza- 
tions (GS and ours) we considered. As a measure of 
the goodness-of-fit we used the confidence level calcu- 
lated from the values of and the numbers of degrees 
of freedom shown in Tables III VII using the CERNLIB 
routine PROB. The agreement of both the GS and our 
parametrizations with both the CMD-2 and SND data is 
on the level that is considered satisfactory in the high en- 
ergy physics (based on /NDF f« 1). Our parametriza- 



tion gives a slightly better confidence level than the GS 
one, especially for the CMD-2 data. 

The disagreement of both parametrizations with the 
KLOE data is disastrous. But the situation radically im- 
proves for our parametrization if wc allow the mass of the 
a;(782) to vary. Adding this free parameter reduces the 
original = 192.6 to = 29.6, what from the statis- 
tics point of view fully justifies this step. The confidence 
level of the hypothesis that our parametrization describes 
correctly the KLOE data has rocketed to 99.62%. This 
is highly non-trivial because the statistical errors of the 
KLOE data arc very small, much smaller than those of 
both sets of data from Novosibirsk. The price that has 
to be paid for this excellent result is lowering the a;(782) 
mass by Am^j = (1.93 ± 0.19) MeV. The physics behind 
this in unclear. We originally suspected that the mea- 
sured values of the mass of the two-pion system might 
contain a systematic shift, at least in the p-w interfer- 
ence region. But the correspondence with the members 
of the KLOE collaboration excluded this possibility. 

Also in the case of the GS parametrization de- 
creased from 367.8 to 202.6 when the u) mass was con- 
sidered an additional free parameter. The required mass 
shift is even higher than in the case of our parametriza- 
tion, namely (2.13±0.35) MeV. Adding the ^-meson con- 
tribution improves the quality of the fit while changing 
the w-mass shift negligibly, compare Tabs. VI and VII. 

On the other hand, the addition of the u) mass among 
free parameters is not required by the CMD-2 and SND 
data. The drop of is less than one both for the GS 
parametrization and ours. Why the KLOE data behaves 
in a different way is unclear. The issue obviously requires 
explanation. We do not have any at the moment. 

Table IX summarizes the results on the p{770) mass 
we have got by fitting various data sets using the GS and 
our parametrization. Its salient feature is that for all 
data sets we considered (CMD-2, SND, KLOE, KLOE 
with the shifted u mass) our parametrization gives a 
smaller p(770) mass than the GS one. The differences 
are large, they range from 7.8 MeV for the KLOE data 
to 10.3 MeV for the CMD-2 data. At the same time the 
errors of the masses themselves are very small. It looks 
like the two columns concerned different quantities. And 
this is probably the core of the problem. The p-meson 
mass is defined in both parametrizations as a parameter 
in the corresponding formula for the propagator. As 
the formulas are different, so must be the definitions and 
the "experimental" values obtained by fitting those two 
parametrizations. 

The cure of this problem is not easy. To get the "uni- 
versal" definition of the p-meson mass it is insufficient 
to declare some formula for the propagator (e.g., that 
based on the GS parametrization) as the canonical one. 
If that preferred formula is not correct, then the different 
processes or even the same process measured in different 
kinematic regions will give different values of the p meson 
mass even if the experiments are ideal. In fact, in order 
to compare the theoretical predictions with the data on 
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TABLE IX: The p meson mass obtained from various data 
sets using the GS and our parametrizations 

GS ours 



CMD-2 

SND 

KLOE 

KLOE" 

KLOE'' 



775.3±1.1 
774.60±0.52 
769.22±0.32 
771.58±0.60 
772.79±0.35 



767.08±0.83 
764.23±0.43 
761.39±0.25 
763.44±0.33 
763.48±0.33 



"Accounting for the energy shift in the p-ui region 
'Accounting for the energy shift in the p-u> region, (/)(1020) in- 
cluded 



TABLE X: The p meson width obtained from various data 
sets using the GS and our parametrizations 



GS 



ours 



CMD-2 

SND 

KLOE 

KLOE" 
KLOE' 



143.2±2.5 
147.3±1.0 
144.71±0.56 

146.39±0.51 
151.53±0.55 



136.1±2.2 
140.29±0.85 
140.07±0.66 

142.56±0.41 
142.78±0.46 



"Fit with the energy shift in the f)u> region 

'Fit with the energy shift in the pu) region, (/> contribution included 



some process, the propagator has to bo convohited with 
other factors ( "vertices" ) that are characteristic for that 
process and depend also on the virtual p-meson momen- 
tum squared. If the propagator is not correct, then the 
cross section or decay width formulas of various processes 
are incorrect, but each in a different way. We illustrate 
this point on a toy example in the Appendix. 

It is generally believed that the only meaningful defini- 
tion of the mass and width of a resonance, which does not 
have asymptotic states, is that in terms of the position 
of the pole on the second Riemann sheet of the propa- 
gator [56]. The pole of the propagator induces a pole in 
the scattering amplitude (at least in perturbative field 
theory). The same propagator enters the amplitudes of 
various processes, so the pole definition should be process 
independent. The pole position does not change if an an- 
alytic function without pole is added to the amplitude to 
describe the background present in experiments. So the 
pole definition should be also background independent. 
In the case of resonances for which a reliable theoretical 
framework exists {W"^ and bosons) is it possible to 
relate the mass and width defined in terms of the pole 
position to the "mass" and "width" that enter as param- 
eters of various fitting formulas [57-59] . The pole-defined 
mass and width are genuine physical quantities (gauge in- 
variant and rcnormalization scheme independent in the 
boson case). 

There have been several attempts to find the position 
of the p^-meson pole on the second sheet of the complex 
s-plane. 



Lang and Mas-Parareda [60] first mapped the cut s- 
plane into a unit disc with data about the P-wavc partial 
amplitude of the tttt scattering falling on the unit circle. 
Then they fitted the data to the product of a pole out- 
side the unit circle (corresponding to the second-sheet 
pole in the s variable) and an analytic function repre- 
sented by its Taylor series. The pole position and the 
coefficients in the Taylor series were parameters of the 
fit. There arc different ways to suppress the importance 
of the higher Taylor series coefficients and make the se- 
ries finite. Unfortunately, this brings some model depen- 
dence. In [60] the Pictarinen [61] smoothness condition 
was utilized. Lang and Mas-Parareda used five different 
data sets and got the mutually compatible values of the 
p mass and width, whereas the values of those quanti- 
ties determined in experimental papers were incompati- 
ble. The pole-defined masses were always lower than the 
masses from the parametrizations used be experimental- 
ists. 

Bohacik and Kiihnelt [62] used the pole-searching 
method developed in [63] on the basis of the statistical 
method for the testing of analyticity [64]. The central 
objects of the latter arc the moments Qn {n natural) cal- 
culated from the data and a smooth function depending 
on the data errors and density. The moments should 
obey the N{Q, 1) distribution if the investigated data de- 
scribe an analytic function. If a pole is present, the mean 
values of the moments are expressed in terms of the pole 
position and residmim. The authors of [62] mapped the 
second sheet on the unit disc, so the assumed pole lay 
inside the unit circle. They either determined the pole 
parameters from the moments calculated directly from 
the data or multiplied the data first by the Blaschke pole- 
killing factor and determined its parameters by requiring 
the calculated Qn momenta be iV(0, 1) distributed. The 
resulting mass and width of the p meson were compatible 
with those of Ref. [60] . This approach suffers from model 
dependence, which crawls in through the construction of 
the error function and the choice how many and which 
moments are taken into consideration [64]. 

Dubnicka and collaborators [23, 29] constructed two 
analytic models incorporating all pion form factor prop- 
erties (including pairs of complex conjugate poles on un- 
physical sheets, corresponding to the p, p' , and p" reso- 
nances) and fitted them to all accessible data. In both 
models, the (pole) mass of the p meson is much smaller 
than the values given by PDG [19]. The authors tried to 
keep the number of free parameters as low as possible. As 
a consequence, the quality of their fits is not very good, 
X^/NDF = 1.81 and 1.58 in [23] and [29], respectively. 

The present work has not been intended as a final word 
in the field of the pion form factor parametrization and 
extracting the resonance parameters. We are even unable 
to offer unambiguous values of nip and Pp. Natural can- 
didates would be those coming from the best fit, which 
is our fit to the KLOE data. But those are hampered by 
the presence of a mysterious energy shift in the p-Lu mass 
region. The main aim of this paper was to issue the warn- 
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TABLE XI: Each group fits its data with its own formula 



Group 


a 


R 


x' 


A 


0.95 


1.096 ±0.016 


11.86 


B 


1.05 


0.866 ±0.013 


10.26 


C 


0.90 


1.369 ± 0.020 


4.88 



ing that the question of the correct parametrizations is 
not simple and requires much attention. The resonance 
parameters obtained by fitting the data are subject of 
large systematic errors of the parametrizations and have 
often been underestimated so far. 
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APPENDIX: A TOY EXAMPLE 



the relative error of 5%. The three generated data sets 
can be accessed at [65] . Each group fitted its data by its 
own formula. The formulas they chose can be written in 
a unified way as 

U = RI", (A.l) 

with a = 0.95, 1.05, and 0.90 for the group A, B, and C, 

respectively. The results of the fits arc shown in Table XL 
The quality of all fits is acceptable according to the high- 
energy physics standards and the resulting parameters R 
have rather small errors. But the Rs of the different 
groups are incompatible. 

Then the groups realized that problem lies in different 
formulas they used and agreed to adopt the formula of 
the B group {a = 1.05) as the common one. They again 
got the good fits and better agreement among their R 



TABLE XII: Each group fits its data with the common for- 
mula (a = 1.05) 



Group 


R x' 


A 


0.928 ± 0.013 10.47 


B 


0.866 ± 0.013 10.26 


C 


0.854 ±0.012 6.10 


TABLE XIII: 


Each group fits its data with the correct formula 


(a = l) 




Group 


R x' 


A 


1.011 ±0.015 4.64 


B 


0.987 ±0.014 8.12 


C 


0.999 ± 0.014 5.18 



Let us assume that three experimental groups A, B, 

and C measured the dependence of the voltage on the 
electric current of the same wire, but in different inter- 
vals. Each interval contained twelve data points. We 
simulated their measurements assuming the validity of 
the Ohm law with R = l^l and the Gaussian distribution 
of experimental values around the Ohm-law values with 



parameters, but some inconsistency, especially between 
groups A and C, remained, see Table XII. 

Only after the groups learnt about the correct formula 
and decided to use common a = 1, they have obtained 
both good fits and compatible results, as it is witnessed 
in Table XIII. 
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